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We describe a simple model of fermions in quasi-one dimension that features interaction- 
induced deconfinement (a phase transition where the effective dimensionality of the sys- 
tem increases as interactions are turned on) and which can be realised using dipolar 
fermions in an optical lattice d. The model provides a relisation of a " soft quantum 
matter" phase diagram of strongly-correlated fermions, featuring meta-nematic, smectic 
and crystalline states, in addition to the normal Fermi liquid. In this paper we review 
the model and discuss in detail the mechanism behind each of these transitions on the 
basis of bosonization and detailed analysis of the RPA susceptibility. 

Keywords: dipolar fermions; optical lattice; phases. 

1. Introduction 

45 years have passed since John Hubbard identified "understanding [...] the balance 
between bandlike and atomic-like behaviour" as the key aim of research in strong 
correlations, and propos ed t he model that bears his name as the simplest embodi- 
ment of that conundrum Yet, in spite of a plethora of systems lying somewhere 
between the Fermi liquid and Mott insulator or Wigner crystal states, the funda- 
mental issue remains unresolved. In fact, even the simple Hubbard model has only 
been solved exactly in one HI and infinte^l dimensions. Given the notable differences 
between the physics of these two extreme cases, one expects very rich behaviour in 
two and three dimensions - as suggested abundantly by experiment. 

Faced with the above difficulty, phenomenological scenarios have been put for- 
ward to deal with the wealth of experimental information on strongly-correlated 
quantum matter. For example, "quantum liquid crystal" phases with partially- 
broken symmetry (or, more generally, incomplete localisation) have been put for- 
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ward as the missing links between the Fermi liquid and crystalline states of the 
fluid of electrons. ® Another, closely related focus of attention have been dimen- 
sional crossovers ^ and, in particular, collective phenomena changing the effective 
dimensionality of the system. (One advantage of focusing on the latter class of 
problems is that the known physics of the one-dimensional Hubbard model can be 
employed as a starting point, and the interaction introduced as a perturbation.) 

In the latter class of phenomena we find the confinement hypothesis. It was orig- 
inally formulated ^ for an array of Luttinger liquids coupled by a transverse single- 
particle hopping amplitude, t±. The hypothesis states that there is a finite, criti- 
cal value of t± below which all coherent motion becomes strictly one-dimensional. 
Indeed recent functional renormalization group calculations ^ show that, for in- 
finitesimally small t±, the ground state with a warped Fermi surface is unstable. 
The instability can be described as an "ironing out" of the Fermi surface, and is 
thus closely related to other interaction- induced Fermi surface shape instabilities, 
notably the Pomeranchuk instability leading into the highest-symmetry of the 
quantum liquid crystalline states: the nematic phased 

Another interesting question to ask is whether the opposite of the confinement 
transition is possible. By this we mean whether the almost flat Fermi surface of a 
system that is quasi-one-dimensional, in the sense that inter-chain hopping is very 
small, can acquire some warping as a result of interactions. A similar phenomenon 
is known to occur in stacks of integer quantum Hall systems In this case, 

the chiral one dimensional Luttinger liquids on the edges of the different layers 
couple together to create a two-dimensional Fermi surface (the chiral Fermi liquid, 
for strong tunnelling and interactions^^ and for weak tunneling ^j) . 

In this article we study this phenomenon in a different theoretical context, 
namely a two-dimensional stack of chains, each of them containing free fcrmions. 
In the absence of interactions the ground state of the system can be described 
as a non-interacting Fermi gas with an almost flat Fermi surface. We introduce an 
inter-chain interaction, and address the stability of the quasi-one-dimensional Fermi 
surface with respect to the perturbation. Interestingly, the inter-chain interaction 
can lead to additional warping, and even the closing of the Fermi surface, which 
thus becomes two-dimensional. In addition, we find that the model also exhibits 
a crystalline state, for sufficiently strong interactions, which is entered through a 
density wave instability. These meta-nematic and crystalline states compete with a 
third, smectic phase corresponding to a different type of density wave instability. 
Thus the model realises some of the phenomenological phase diagram mentioned 
above 

In the following sections we introduce the model, and discuss how it could be 
realised using cold atoms We will then describe in detail the mechanism by which 
the different phase transitions in the model come about, using a detailed analysis 
of the RPA susceptibility. We also offer a discussion of some aspects of the physics 
of the system using arguments based on bosonization. 
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Fig. 1. (a) Our Hamiltonian describes hopping mostly along a set of chains (with very little 
perpendicular hopping) plus a nearest-neighbour, inter-chain interaction, (b) A two-dimensional, 
strongly-anisotropic optical lattice with dipolar interactions can provide an approximate experi- 
mental realistion (see text). 

2. Model 

The unperturbed Hamiltonian is 

^ho P = -*n £ Ia,i+i - ^XX^+m + H - c - - mXXa*. oo 

l : i l : i l,i 

It describes a set of chains kept at the same chemical potential, fi. c| • creates 
a fermion on the site of the l™ 1 chain. The particles can hop easily along the 
chains, with amplitude tii , but not so much between the chains, for which the hopping 
amplitude is t± <C tii. We assume the spin degree of freedom to be frozen (spinless, 
or fully polarised, fermions -see below). 

The above Hamiltonian is almost diagonal in the chain index, I, and therefore 
has a nearly- flat, quasi- ID Fermi surface. 

The perturbation is an interaction between particles sitting on adjacent sites of 
different chains: 

Ant = ( 2 ) 

l,i 

The full Hamiltonian, 

ff = #hop+- H int> ( 3 ) 
is represented schematically in Fig. [T] (a) . It describes hopping along the chains (the 
solid lines in the figure), with only small inter-chain hopping, plus an inter-chain 
interaction (represented by the dashed lines). 

3. Experimental Realisation 

A ne w av enue into strong correlations has been opened recently in ultracold atomic 
gases Indeed the Mott transition has been observ ed in v ery precise experimental 
relisations of the Hubbard model in a harmonic trap H41H1 i n principle, a sequence 
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of experiments on systems with different trap profiles could be used to deduce the 
phase diagram of the Hubbard model in the absence of the trapQ^. Moreover, even 
simpler models can be realised to explore the physics of dimensional crossovers and 
phenomenological scenarios of strong correlations. Interestingly, our model can be 
realised, quite precisely, in an o ptical lattice setup loaded with quantum-degenerate 
dipolar fermionic molecules yjJ or fermionic atoms with large magnetic dipoles such 

The Hamiltonian ([3|) is theoretically convenient because the unperturbed part 
features no interactions and can therefore be described in terms of free fermions. 
Interactions are introduced perturbatively. Moreover the interactions do not change 
the nature of each individual chain, but take place between distinct chains exclu- 
sively. In this section we describe a possible realisation of this physics. 

The atoms (or, equivalentely, molecules - though for concreteness we will assume 
atoms in what follows) are trapped in a two-dimensional optical lattice by two pairs 
of counter- propagating laser beams: see Fig.[T](b). In the figure, the solid and dashed 
lines represent their respective wavefronts. If the lasers are sufficiently intense, the 
lattice can be described in the tight-binding limit, with only one orbital per site. 
Moreover by making one of the pairs of lasers much more intense than the other 
we can create "chains" along which hopping can take place, but such that hopping 
between different chains can be very small. These chains are represented by the 
solid lines in Fig.[T](b). The lattice sites are where the solid and dashed lines cross. 

To suppress the interaction between atoms that are in the same chain, we exploit 
the dependence of the dipole-dipole interaction, 

t ,/-~n ,9 1 — 3cos 2 # 

( )= |R| 3 ' () 

on the angle 9 between the vector giving the relative positions of the two dipoles, R, 
and the external field, H. This interaction is identically zero for 9 — arccos (l/ v3j ~ 
54.736°. Thus aligning the chains at this angle (as in the figure) we ensure that there 
are no interactions between atoms in the same chain (on-site interaction, including 
the additional one due to the Van der Waals forces between the atoms, is already 
suppressed by Pauli's exclusion principle, for fully polarised fermions). 

Finally, by arranging the chains so that each site is on the horizontal from the 
closest sites on the two adjacent chains we ensure that the interaction with those two 
sites is maximally repulsive [V (R) = d 2 / |R| 3 , corresponding to 9 = tt/2]. Longer- 
ranged inter-chain interactions with other sites on the two adjacent chains can be 
made comparatively weaker by making the lattice constant very long along the 
longitudinal direction. This is aided by the relatively rapid fall of the interaction 
potential with distance, cx 1/|R| . To achieve this in combination with keeping 
hopping along the dotted lines to a minimum requires that the laser creating the 
solid line wave fronts be very intense. This ensures a very high potential barrier 
for the atoms to hop between chains. It is also desirable to orient the plane of the 
lattice so as to ensure that all inter-chain interactions are repulsive. The details of 
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how the precise realisation of our Hamiltonian can be achieved are given in Ref. 
4. Meta-nematic phase transition 

In order to establish which ground states may occur in this system we start by 
evaluating the stability of the Fermi surface shape. The interaction in Hamiltonian 
^ wants to avoid having two particles on the same site on neighboring chains. 
Physically, this is very similar to interchain hopping - if only one electron is present, 
it can lower its kinetic energy by hopping. We may therefore expect the interaction 
to renormalize the effective interchain hopping to be larger than its bare value. 
We investigate this possibility by using a res tricted Hartree-Fock mean field theory 
similar to those used to study Pomeranchukl^2^ and topological^ Fermi surface 
shape instabilities. 

We use as a trial ground state a Slater determinant of plane waves, 

l*)=n[( 1 -^) + A r k 4]|0), (5) 

k 

determining the occupation numbers = 0, 1 by requiring that the momen- 
tum distribution minimizes H ^). The momentum distribution corre- 
sponds to a non-interacting Fermi gas with a renormalized dispersion relation 
e£ = — 2t\\ cos (fc||) — 2t* L cos(k±) — fJL* . The structure of the interaction H[ n i is 
such that only the perpendicular hopping is changed. It is given by the following 
self-consistency equation: 

t* ± = t± + ^^cos(fc ± )7V k , (6) 

k 

where = 9(— e£) is the Heaviside function. In what follows all references to the 
chemical potential will be to its renormalized value so we will omit the * for that 
quantity. 

Numerical solutions to this equation are shown in Figs. [2] and [3] As either 
the bare inter-chain hopping t± or the interaction V is increased, the renormal- 
ized hopping, t* ± , initially increases linearly but then has two bifurcation points, 
between which lies a first-order jump to a higher value. The order parameter of 
this phase transition is the amount of delocalisation in the perpendicular direction, 
ip = (cjjCj^+i + H.c). We refer to the jump of ip as we vary t± as a meta-nematic 
transition in analogy with meta-magnetism (where the magnetisation jumps under 
an applied magnetic field). 

Focusing on Fig. O we now discuss some features of this transition. For fixed /i, 
the transition always occurs at the same value of the renormalized hopping tj_ , and 
in fact corresponds exactly to the point where the Fermi surface switches from open 
to closed, something we refer to as dimensional crossover from quasi-lD to quasi-2D. 
We see also as we take the Fermi level from half-filling /i = to nearer the bottom 
of the band /i = — 1.5tii the transitions generally moving down to a lower energy. In 
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some of the lines, there are no transitions - which is when the bare perpendicular 
hopping t± is already large enough for the bare Fermi surface to be closed - and the 
renormalization occurs in a continuous manner. Another point to note is that in 
some of the lines, there are two separate first order transitions - the second in fact 
corresponds to the value of t± where the Fermi-surface becomes open again, but 
in the other direction. This gives strong hints that the meta-nematic transition is 
indeed a dimensional crossover phenomena - however this second transition requires 
very large renormalization of tj_ well outside any regime where Hartree-Fock is valid, 
and furthermore one would expect the electrons to crystallize before that point (see 
the next section) so we will not say anything further about the second transition. 

As V — > 0, the meta-nematic transition becomes more and more weakly first 
order and requires a larger bare value of t±. At strictly V = there is no longer 
a first-order phase transition, but the phenomenon survives at t±_ = t\\ + /i/2 as 
a 'two-and-half order Lifshitz transition^ - see Fig. Ufa). On the other hand the 
transition is first order for any V > 0. 

The meta-nematic transition results from enhanced scattering when the poten- 
tial reaches the singularities at the edges of the ID bands. This is a density of states 
effect - see Fig. H[b) and hence we expect it to be robust to quantum fluctuations 
present for large values of V/t\\ and not taken into account by our mean field theory. 



5. Crystallisation 

The meta-nematic transition is not the only one possible in the system described 
by Eq. j3]). In fact, for t± = 0, the Fermi surface is perfectly nested, leading to the 
possibility of a density wave instability at low temperatures - indeed, the strong 
coupling phase of Hamiltonian [3] is a crystalline state. Another way of thinking 
of this is in terms of the interchain backscattering- - . Although in the cold atom 
set-up, the fermions are not charged, we will still refer to this as a charge density 
wave (CDW) instability in order to draw parallels with previous work on electronic 
models. 

We probe the potential CDW instability by examining the Fourier transform of 
the dynamic susceptibility, 

X(M)=i<¥|Tp(M)pt(k,0)|*). (7) 

Here p(k) = ^ q c q c q _k is the Fourier transform of the local occupation number, 
and p(k, t) is its Heisenberg representation. 

The 'noninteracting' susceptibility (Lindhard function) is given by 

A (k,w)= / — — 2 — « , ♦ — ■ (8) 
J (2tt) 2 uj - e q + e q+k 

We note that this is written in terms of the renormalised dispersion relation e£ 
(technically the 'bubble' should be dressed). We also point out that because the 
major contribution to the kinetic energy is perpendicular to the direction of the 
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Fig. 2. The renormalized hopping amplitude tj_ as a function of interaction strength V for various 
different values of chemical potential fi and bare hopping tj_ . The values of V where tj_ is not single- 
valued which ultimately corresponds to a jump in as V is increased, mark the meta-nematic 
transition. A full discussion of these graphs is given in the main text. 




Fig. 3. A different way of viewing the data: The renormalized hopping amplitude tj_ as a function 
of the bare hopping t± for various values of interaction strength V and chemical potential fi. 
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Fig. 4. (a) A plot of the renormalized hopping t*, as a function of interaction strength V for 
different chemical potentials fi approaching the bottom of the band when the bare hopping t± = 0. 
This figure clearly shows both the critical value of V and the size of the first order jump decreasing 
towards zero as the bottom of the band is reached - i.e. the Lifshitz transition, (b) The density of 
states of the anisotropic 2D lattice for various values of the ratio between parallel and perpendicular 
hopping. The van-Hove singularity in the density of states occurs at exactly the energy when the 
contour of constant energy just touches the edge of the Brioullin zone - i.e. when the Fermi surface 
switches from open to closed. Enhanced scattering due to this singularity is the origin of the 
meta-nematic transition. The small tails at the band edges are a numerical artefact. The energy is 
measured in units of t||. 



interaction, any vertex corrections to the bubble are small and may be neglected in 
lowest order. 

The interaction may then be treated within the random phase approximation 
(RPA) 

yf] \ _ A Q (k,^) /Q\ 

X[K l + 2Vcos(fc ± )Ayk, W )- (9) 

An instability at wave vector k takes place if the static component of the suscepti- 
bility diverges, X(k\\,k±,u) = 0) — > oo. 

For t* ± <C ill the Lindhard function is strongly peaked at (2fejr,7r), and in fact 
logarithmically divergent at this wavevector when t* ± = - see Fig. ■ This is 
due to the strong nesting of the quasi one-dimensional Fermi surface, and it makes 
the system unstable to a CDW of that periodicity at a critical coupling V given by 
the following Stoner criterion: 

1 = 2V c X (2k F ,Tr,uj = 0), (10) 

which clearly gives V c = in the perfectly nested Fermi-surface when t± = 0. More 
generally one has to evaluate Xo(2kp,ir,ui = 0) to obtain V via Eq. (flQ|) . We first 
plot the results in FigE^b). The physical picture is then as follows: begin at V = 
where the system is a quasi-lD Fermi gas. Start increasing V, and the interchain 
hopping will renormalize to its value t* ± . As V is increased further, eventually the 
CDW instability Eq. (|10| will be satisfied, with Xq evaluated using the renormalized 
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hopping fj_. At this point, the system becomes crystalline, and will remain crys- 
talline as V is further increased arbitrarily. If this occurs after the meta-nematic 
transition, then both phases will be seen. If the crystallization happens first (as is 
the case when t± = 0), then the meta-nematic transition will not be seen. This 
allows the full phase diagram of the Hamiltonian ([3]) to be plotted in Fig. [6] 




Fig. 5. (a) The Lindhard function, Eq.[8]as a function of (fen , k±) along the path (tt, 0) — > (0, 0) — > 
(0, tt) — » (tt, tt). For clarity, the axis has been shifted vertically for the three different values of t± 
plotted. The singularity at (2k p , tt) is clearly seen in the case when t± = 0, reducing to a sharp 
peak when tj_ is finite, (b) The critical V c from Eq. 1101 plotted against the bare hopping t± for 
various values of fi. The humps are at the points when the Fermi surface switches from being open 
to closed - which results in much less nesting. 



6. Smectic phase 

So far, we have been using the nearest neighbor approximation for the interaction 
- Hamiltonian [3] In this case, the Fourier transform of the interaction potential 
V(k) is independent of fcn, hence the dominant CDW instability is always at the 
peak of the Lindhard function, i.e. fc|| = 2kp. However, we now consider the full 
structure of the dipole interactional In particular, when the ratio of lattice spacings 
a — a\\/a± is not large, T^(k) acquires a large dependence on For small values of 
the parameter a, the largest negative value of V(k) is at k = (0,7r). Although the 
Lindhard function itself is small at this wavevector (Fig. [5]), the product V r (k)X (k) 
may not be, and can become larger than the product VX Q at (2kp,ir). Hence so 
long as t* ± (i.e. Xo(k\\ = 2kp) is finite), then there is a level of anisotropy of 
the lattice where the leading instability is at (0,7r) and not (2/cp,7r). 

The (0, 7r) instability is still a form of CDW - however as it breaks lattice sym- 
metry in one direction only, it has smectic order. Fig. [7] shows which of the two 
instabilities takes place first as the strength of the interaction is increased. The 
smectic order is favored when the interaction between neighboring chains is such 
that the fermions can lower their energy by crowding every other chain, paying a 
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Fig. 6. Phase diagram of the Hamiltonian of Eq. Q for fi = — 1.9tii (a) and — l.Stii (b). The 
circles track the two bifurcation of the first-order meta-nematic quantum phase transition between 
the quasi-one dimensional (q-lD) and two dimensional (2D) phases. The solid line marks the line 
of quantum critical points separating these two Fermi liquid states from the crystalline state (C). 




Fig. 7. The phase diagram that shows whether the dominant instability is towards the crystalline 
(C) or smectic phases (S). The axes are the lattice anisotropy a^/a± and the relative strength of 
the perpendicular hopping tj_, for n* = ±1.9tii (a) and ±1.5tii (b). 



penalty in kinetic energy but lowering the interaction energy. Note that the strong- 
coupling limit ground state is always a density wave. 
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7. Bosonization Approach 

In the particular case t± = 0, the 'bare' dispersion relation is strictly one- 
dimensional, allowing us, in principle, to employ the bosonization technique. In fact, 
the bosonized Hamiltonian features backscattering terms that are responsible for 
the checkerboard crystallization which, in this limit, takes place at arbitrarily small 
coupling Thus strictly speaking a bosonized description of the fluid of fermions 
in this model is never valid. Nevertheless it is illustrative to ignore the backscatter- 
ing terms and see what happens. The Hamiltonian can then be diagonalized exactly 
to compute the holon dispersion relation in the abscence of crystallization. 

To bosonize the interacting part of the Hamiltonian first we do a Fourier trans- 
formation of all the creation/annilation operators and then we manipulate the order 
of the operators further to produce a form where only the fermionic density appears. 
Note that we will need to work in two-dimensional space since the interactions cou- 
ple the chains. At the end the interacting part of the Hamiltonian becomes: 

Hint = - 2V E cos(q x )p(q = 0) + 2V £ cos( fe )p(q)/S(-q) (11) 

q,9j;>0 q,(J x >0 

where V is the potential and the direction x is the one perpendicular to the chains. 
The first term then shifts the chemical potential since p(0) is the density operator 
at q = (in other words it counts the number of fermions). 

The most interesting four-fermion term then, provides in the language of Ref . ^31 
: g' 2 = .94 = 2V cos(q x ) / (2ttvf) where vf is the Fermi velocity. Then the dispersion 
relation becomes u p (q}\q v \, with: 



Wp(q) = Vf\ 



l 2Vcos{q x ) 

■KVp 



-I 2 


~2Vcos(q x )~ 




TTVp 



(12) 



If the expression under the square-root becomes negative then the energy of 
low-lying excitations becomes unphysically imaginary, signalling the breakdown of 
bosonization. The condition for this to take place is: vf > —^Vcos(q x )/ir which, 
if it holds for every q. x , demands ir/Avp > V. Thus we find that the holon velocity 
vanishes at a critical coupling which corresponds to the onset of the smectic phase 
and coincides exactly with a divergence of the RPA susceptibility at (tv,Q) at the 
same critical coupling V cr it- Beyond Verity the bosonized theory has broken down 
and it is no longer useful. Moreover we stress that even below this value the use of 
bosonization is somewhat artificial, since we need to cross out the backscattering 
terms that, in reality, will lead to the checkerboard crystallisation before any of 
the physics that we have discussed on the basis of that approach take place (as our 
RPA approach shows). However it affords a different perspective on the nature of the 
transition into the smectic state. In particular, it shows that the smectic instability 
can be interpreted in terms of plasmons going soft as we approach it. Moreover, 
it lends further credibility to our RPA approach, which can be applied over the 
whole phase diagram and used to find both the instabilities into the smectic and 
crystalline states, as the critical coupling V cr it obtained with bosonization for the 
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smectic phase coincides with that given by bosonization, in the appropriate limit 
t± — > (although, we stress once again, in that limit the crystallisation takes place 
first). 

8. Conclusions 

With the help of optical lattices and cold atoms we can realise correlated fermion 
physics. In particular interactions between dipolar atoms or molecules can lead to 
quantum soft-matter phases, hypothesized in strongly interacting quantum systems. 
Meta-nematic, smectic, and crystalline phases are predicted to occur for fermions. 
For bosons the expected phases are either a superfluid phase or a CDW. 

The tunability of the experimental parameters in optical lattices leads to many 
possibilities for experiments with polarised fermions. For example, by orienting the 
field along the chain direction, we could realise a model with attractive intra-chain 
interactions and repulsive inter-chain interactions, presumably leading to supercon- 
ducting stripes. Using the same set-up discussed here, but with the chains perpendic- 
ular to the direction with no interactions, we can have one-dimensional Luttinger 
liquids and then, by adding transverse hopping one can study the fundamental 
confinement-deconfinement transition and orthogonality catastrophe. 

Another possible extension of the present study is the three-dimensional optical 
lattices as well as the regime of finite temperatures. 
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